Composite-fermionization of bosons in rapidly rotating atomic traps 
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The non-perturbative effect of interaction can sometimes make interacting bosons behave as 
though they were free fermions. The system of neutral bosons in a rapidly rotating atomic trap 
is equivalent to charged bosons coupled to a magnetic field, which has opened up the possibility 
of fractional quantum Hall effect for bosons interacting with a short range interaction. Motivated 
by the composite fermion theory of the fractional Hall effect of electrons, we test the idea that the 
interacting bosons map into non-interacting spinless fermions carrying one vortex each, by com- 
paring wave functions incorporating this physics with exact wave functions available for systems 
containing up to 12 bosons. We study here the analogy between interacting bosons at filling factors 
v = n/(n + l) with non-interacting fermions at v* = n for the ground state as well as the low-energy 
excited states and find that it provides a good account of the behavior for small n, but interactions 
between fermions become increasingly important with n. At v = 1, which is obtained in the limit 
n — > oo, the fermionization appears to overcompensate for the repulsive interaction between bosons, 
producing an attractive interactions between fermions, as evidenced by a pairing of fermions here. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sation (BEC) of atomic gases^ 2 - has generated a rich va- 
riety of phenomena. In particular, it has allowed the pos- 
sibility of testing the remarkable concept of "statistical 
transmutation," namely the idea that interacting bosons 
may sometimes behave like spinless fermions. For con- 
tact interactions, it may seem rather sensible for bosons 
to emulate fermions, to the extent allowed by symme- 
try requirements, because the Pauli principle itself fully 
takes care of the repulsion. Of course, a conceptual un- 
derstanding of how this precisely happens, what it means, 
and how bosons can behave like fermions while satisfying 
the constraints of bosonic exchange symmetry requires 
a detailed theory. The tendency for fermionization has 
been appreciated for quite some time for bosons in one 
dimension^i4 Girardeau 3 , showed that for an infinitely 
strong delta function repulsion, the bosonic ground state 
wave function \t B is related to the Slater determinant 
ground state wave function \E' F for spinless fermions in 
one dimension as: 
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(1) 



The problem was solved exactly for an arbitrary strength 
of the interaction by Lieb and Liniger the fermions de- 
scription is a useful starting point in the strong-coupling 
limit, when the interaction strength is large compared 
to the Fermi energy. Recent experiments 5 ^ are in excel- 
lent agreement with the Lieb-Liniger theory in the entire 
range of interaction strength, which can be varied in an 
optically confined one dimensional boson system by con- 
trolling the density and the confinement strength. 

This work is concerned with the possibility of an emer- 
gence of fermion-like structures in a bosonic system in 
two dimensions, under conditions appropriate for a frac- 
tional quantum Hall effect (FQHE) of bosons. The famil- 



iar FQHE occurs when charged electrons are confined to 
two dimensions and exposed to a strong magnetic fieldi 
There is no realizable system of charged bosons where 
FQHE can be studied. However, a system of neutral 
atoms in a rotating trap is mathematically equivalent to 
a system of charged bosons in a magnetic field, which, 
with confinement to two dimensions, should create, for 
sufficiently rapid rotation, a FQHE state of bosons. BEC 
systems confined to two dimensions have been created^ 
and their properties have been studied under rotation, 
although the FQHE conditions have so far not been 
achieved. Rotation of a BEC produces vortices in the 
condensateiSiiSiii As the rotation frequency is increased, 
the BEC state is destroyed and, eventually, the FQHE 
state may be achieved (the latter has no off-diagonal long 
range order). These advances have motivated a number 
of studies of the FQHE of bosons interacting via a short- 
range interaction, 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 

We will assume below that the Landau level (LL) spac- 
ing for bosons is sufficiently large that it is a good ap- 
proximation to restrict the bosons to the lowest Landau 
level. The bosonic system is then always in the strong 
coupling limit, because the nature of the state is com- 
pletely determined by the interaction. In fact, the solu- 
tions are independent of the strength of the interaction, 
which merely sets the energy scale. It is natural to ap- 
peal to the fractional quantum Hall effect of electrons for 
guidance. Laughlin's wave functional can be generalized 
for the ground state at the bosonic filling v = 1/2: 
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where Zj = Xj — iyj denotes the position of the jth 
boson on the two-dimensional plane, and the magnetic 
length has been set to unity. More generally, the un- 
derstanding of the electronic FQHE is based on the for- 
mation of quasiparticles known as composite fermions 
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(CFs); specifically, the sequence of fractional filling fac- 
tors v e = n/(2n + 1) (y e refers to the electronic filling 
factor) is understood as the integral sequence v* — n of 
composite fermions i^SjS 3 . Application of completely anal- 
ogous ideas raises the possibility that interacting bosons 
at v = n/(n+l) may behave like free fcrmions at v* = n. 
Jain's wave functions^ can be generalized to 

=V LLL \{{z J -z k )^ (3) 

j<k 

where is the wave function for non-interacting 

fermions (at the effective filling factor), and Vlll 
projects the wave function into the lowest Landau level. 
An explicit mapping between interacting bosons and non- 
interacting fermions should be noted. Eq. produces 
wave functions for the ground and excited states at ar- 
bitrary filling in the range 1 > v > 1/2. This paper 
examines their accuracy by comparison with exact wave 
functions. If valid, a simplification of the problem is 
achieved through a mapping of a non-trivial interact- 
ing boson problem into a more amenable non-interacting 
fcrmion problem, and many essential properties of bosons 
in rapidly rotating traps should find an explanation in 
terms of almost free particles. 

In addition, at filling factor v = 1, we will consider 
Moore and Read's Pfaffian wave function 2 ^, given by 




PfjMjfc} is the Pfaffian of an antisymmetric matrix M 
with elements Mjk, defined as (up to an overall constant) 

Pf{Af ife } = ^2 sgn(o-)M CT (i) .(2)M .(3) CT ( 4 ) . . . M^jv-i)^), 

(5) 

where the sum is over all permutations <r, sgn(er) is +1 
or —1 depending on whether the permutation is even or 
odd, and N is an even integer. The Pfaffian has the same 
form as the projection of the real space Bardeen-Cooper- 
Schrieffer wave function into a fixed number of particles 
N, and therefore represents a paired state of fermions, as 
noted by Greiter, Wen and WilczekiS The fermion pair- 
ing manifests through an incompressible state of bosons. 

The mapping into fermions for the bosonic FQHE 
problem is conceptually distinct from that applicable in 
one dimension (Eq. (JTJ). The modulus of the fermion 
wave function is a manifestly bad approximation for the 
former, because such a wave function has substantial 
mixing with higher Landau levels, and therefore a very 
high kinetic energy. 

Much work has already been done toward testing the 
composite fermion theory for interacting bosons in a mag- 
netic field. Many studies take bosons to be in a plane, 
confined to a disk by a parabolic confinement; these 
are analogous to the CF theory of electrons confined 
to a parabolic quantum dotiSS*2i Viefers, Hansson and 
Reimannji^ Cooper and Wilkin 13 , Wilkin, and Gunni^, 



and Manninen et alJl have found high overlaps between 
the exact solutions and Jain's wave functions for up to 
N = 10 particles at the "magic" angular momenta of the 
yrast spectrum; further, they also found that the state at 
v = 1 is well described by Moore-Read's wave function. 



While a parabolic potential appears naturally for op- 
tically confined bosonic systems, the strength of confine- 
ment can be varied, and it may be useful to consider the 
situation without confinement. For a large number of 
bosons, it is natural, in the simplest approximation, to 
neglect the effect of boundaries and concentrate on the 
bulk properties. That is most conveniently accomplished 
in theory by studying bosons in the spherical geometry, 28 
in which the bosons move on the two-dimensional surface 
of the sphere, with a radial magnetic field produced by a 
magnetic monopolc at the center. Exact diagonalization 
studies have been carried out in the spherical geometry. 
Regnault and Jolicoeur— have shown that the ground 
state at v = n/{n + 1) is incompressible, consistent with 
the analogy to filled LL state at v* = n. Their results 
also show evidence of incompressibility at v = 1. These 
authors 2 ^ and Nakajima and Ueda— have studied the 
excitation spectrum at v = 1/2. Xie et al? . had earlier 
studied charged bosons in the lowest LL; we will briefly 
consider bosons with long range Coulomb interaction at 
the end, but our focus will be on bosons interacting with 
a short-range interaction, as appropriate for the atomic 
system. None of these studies, however, has carried out 
a microscopic comparison of the exact eigenstates with 
the wave functions of Eq. Another possible geome- 
try without boundaries is the toroidal geometry, which 
has been employed by Cooper, Wilkin and Gunni& in the 
context of bosons in rotating traps. 



We will consider below the spherical geometry and re- 
port on detailed and quantitative tests of the validity 
of the correspondence between interacting bosons in the 
FQHE regime and free fermions in the integral quan- 
tum Hall regime, which makes definite predictions for the 
quantum numbers of the low-energy states of the inter- 
acting boson system, their energies and their cigenfunc- 
tions. Various trial wave functions will be compared with 
the exact eigenstates and the predicted energies with the 
exact eigenenergies. It is well known that Laughlin's 
wave functionjSi which is also a special case of Eq. (0, 
is the exact solution for the ground state at v = 1/2 
for bosons in the lowest Landau level interacting with a 
short range interaction. However, that by itself does not 
imply a correspondence between interacting bosons and 
free fcrmions; for that purpose it is necessary to verify 
the correspondence of Eq. (J3J) for the ground states and 
excitations in a broader range of filling factors. We will 
test it for the ground state and excitations at v = 1/2, 
v = 2/3 and v = 3/4. 
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II. THE HAMILTONIAN 

We consider a system of N bosons with mass m in a 
harmonic trap that is rotating with frequency uj. In the 
rotating frame of reference, the system is described by 
the Hamiltonia n? 1 ^ 2 




<i — mojz x Ti) 2 



+ |[(^- W 2 )(^+ y4 2 )+^]} 

JV 

+ Y,V(r i -r j ), (6) 

i<j 

where ui r and u> z are the radial and axial trap frequencies 
respectively. Vectors = (xi,yi,Zi) represent particle 
positions. In an ultra-cold dilute Bose gas, the scattering 
between particles is dominated by the s-wave scattering 
process. It is then an excellent approximation to describe 
the interaction by a delta function^ 

n ,) = ^)(,), (7) 
m 

where a s is the s-wave scattering length, assumed to be 
positive in this work. When uj and ui r are identical, the 
Hamiltonian resembles that of particles with charge e in 
a magnetic field B = (2muj/e)z. An effective magnetic 
length is defined as £ = y/h/(2muj). The effective cy- 
clotron frequency is defined as uj c = eB/m — 2uj. If the 
axial trap is strong enough such that the wave function 
along the z direction is the ground state of the harmonic 
potential in the z axis, the system enters a two dimen- 
sional (2D) regime where the potential felt by particles 
is written as2S 

V(r)=g5^(r), (8) 

with g = h 2 a s ^/&TT / (ml z ) , where l z — \J h/(muj z ). The 
energy scale in the 2D regime is set by the effective cou- 
pling constant g. We will assume below that the interac- 
tion strength is sufficiently small compared to the Lan- 
dau level spacing that LL mixing is negligible. From our 
experience with electronic FQHE, we know that a mod- 
est amount of LL mixing does not significantly alter the 
results. 



III. COMPOSITE FERMION THEORY 

For bosons in the lowest Landau level, there are three 
situations, (i) For v < 1/2, there are many linearly inde- 
pendent wave functions that vanish upon coincidence of 
bosons, producing an enormous ground state degeneracy, 
(ii) For v = 1/2 there is a single wave function that has 
zero energy for the delta function interaction, giving a 
non-degenerate ground state here. It remains the ground 



state for arbitrarily high coupling g, and may be consid- 
ered to be the analog of the Girardeau wave function of 
the one dimensional problem, (iii) For the excitations at 
v = 1/2, or for any eigenfunctions at v > 1/2, there are 
no wave functions in the lowest Landau level that vanish 
when two particles coincide. While no exact results are 
available here, analogy to fermions gives plausible wave 
functions that we now describe. 

A. u = l/2 

For v = 1/2, Laughlin's wave function for the ground 
state is given by (in the spherical geometry) 

*f/ 2 = *i, (9) 

where 

$i = _Q(ttjUfc - VjUk), (10) 

j<k 

is the wave function of the lowest filled Landau level. 

(uj,Vj) = (cos(6' J /2)e* /2 , sin(6» J /2) e -* /2 ) (11) 

are the spinor coordinates describing the position of a 
particle on the surface of a sphere. It is the exact ground 
state for bosons at v = 1/2 interacting with a delta func- 
tion interaction, which can be seen straightforwardly by 
noting that <I> 2 is the only wave function at v = 1/2 
that is confined to the lowest LL and has zero interac- 
tion energy for the delta function interaction. The wave 
functions $ 2p with p > 2 are not relevant for the short 
range interaction, as these are degenerate with a large 
number of other states. 

B. ^ > 1/2 

For electrons in the lowest LL, the CF theory 2 ^ hy- 
pothesizes that strongly interacting electrons map into 
weakly interacting fermions of a new kind, called com- 
posite fermions. The composite fermions experience an 
effective magnetic field given by B* = B — mpfio, where 
B is the external magnetic field, cf>o = hc/e, and m is 
an even integer. Equivalently, the filling factor of com- 
posite fermions, v* , is related to the electron filling fac- 
tor, v, by v = v* /{mv* + 1). This is interpreted in 
terms of electrons having captured an even number (m) 
of flux quanta of the external magnetic field to become 
composite fermions, which no longer see the magnetic 
flux that they have assimilated into themselves. This 
physics suggests the wave functions "ff = Vlll&T^v* > 
where $> v * is the Slater determinant wave function for 
non-interacting electrons at i/*, $i is the wave function 
of one filled Landau level, and Vlll projects the wave 
function into the lowest Landau level, as appropriate for 
very large magnetic fields. These wave function explic- 
itly relate the eigenfunctions of interacting electrons at 
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v to those of non- interacting fermions at v* , and have 
been tested both in the spherical geometry^ and the 
disk geometryi2SiSL24 

The considerations in the preceding paragraph are 
readily generalized to bosons by taking m to be an odd 
integer. We specialize to m — 1 (other odd integer val- 
ues not being relevant to the problem of our interest) and 
filling factors 



which correspond to v* = n of fermions. The wave func- 
tion at v is now given by 

=V L LL^n- (13) 

which is the spherical analog of Eq. ©. The Jastrow 
factor $i now attaches a single vortex to each fermion in 

These two equations define the mapping between in- 
teracting bosons and non-interacting fermions in micro- 
scopic detail. The first equation has implications about 
the structure of the low-energy eigenstates of the inter- 
acting boson system, whereas the last gives trial wave 
functions for the eigenstates, and also the eigenenergies. 
There are two ways to physically think about the above 
equations, (i) Bosons have captured an odd number of 
vortices each to convert into a composite fermion^ (ii) 
The bosons are represented as bound states of fermions 
and an odd number of vortices. The tests below, of 
course, are independent of the interpretation. 

We note that the ground state and excitations of inter- 
acting bosons at v = n/(n-\- 1) are images of the ground 
state and excitations of fermions at v* = n according 
to Eq. (|13|l . The wave function for the ground state at 
v = n/(n + 1) is given by Eq. (|15|) with <&„ taken as 
the wave function of the ground state at v* = n, i.e., 
the n filled Landau level state. The wave function for 
the excited state is similarly related to the lowest energy 
particle-hole excitation, i.e. an exciton at v* = n. The 
eigenstates of the spherical geometry are labeled by the 
total orbital angular momentum, L. The ground state 
has L — 0, which implies uniform density. It has no 
adjustable parameters, given that the wave function of 
n filled Landau levels is unique. The wave function for 
the exciton for any given L is also determined completely 
by group theory, and therefore is free of any adjustable 
parameters. 

A subtle feature of the composite fermion theory ought 
to be noted. The ground state wave function at v = 1/2 
{v* = 1) is given by \frf/ 2 = ( n0 lowest-Landau level 
projection is required here, because the wave function is 
already in the lowest Landau level). It manifestly elim- 
inates spatial coincidence of particles, and thus has zero 
interaction energy for the contact interaction potential. 
As mentioned earlier, no such wave functions can be writ- 
ten, even in principle, for the excited states at v = 1/2 
or for any states at v > 1/2. The CF theory circumvents 
this problem by first neglecting the lowest LL constraint 



to write wave functions ($i$„) in which bosons do not 
occupy the same spatial position, and then projecting 
them into the lowest Landau level, hoping that this would 
capture the actual correlations within the lowest LL. The 
wave functions *f? B are in general much more complicated 
than Laughlin's wave function at v — 1/2. Their validity 
is far from obvious, and their confirmation would provide 
a non-trivial evidence for the composite-fermionization of 
the bosonic system. 

C. v = 1 

We will also be interested in the nature of the state in 
the limit of n — > oo, i.e. at v = 1. Let us recall what 
happens for electrons in this limit, which corresponds to 
v e = 1/2 for electrons. If the residual interactions be- 
tween composite fermions are negligible, a Fermi sea of 
composite fermions is obtained here (the state with an 
infinite number of filled Landau levels is another repre- 
sentation of a Fermi sea), as proposed by Halperin, Lee, 
and ReadiSi That provides a good description of the com- 
pressible state at v e = 1/2, also explaining why there is 
no FQHE here. 3 X However, in the second Landau level, 
electrons form an incompressible state when the Landau 
level is half full (which corresponds to a total filling of 
v e = 5/2), which appears to be best described by Moore- 
Read's wave function. This implies that the mapping 
into noninteracting composite fermions is no longer valid, 
and one must consider the residual interactions between 
them, which presumably cause a pairing instability of the 
CF Fermi sea^ 

If bosons behaved like non-interacting fermions in the 
limit of n — ► oo, the system at v = 1 would be analo- 
gous to the Halperin-Lee-Read Fermi sea. On the other 
hand, if the bosons map into interacting fermions, Moore- 
Read's wave function becomes a plausible candidate. In 
the spherical geometry, it is given by 




IV. CALCULATION 

We will study the wave function in Eq. (|13fl for the 
ground states and excitations at v = 1/2, v = 2/3 and 
v = 3/4. For technical convenience, we will define the 
lowest LL projection as follows: 

*f/(n+l) = ^LLL^n (15) 

The philosophy of "lowest LL projection" is to obtain 
from the unprojected wave function "I ) i < I ) n , which has 
some amplitude in higher Landau levels, a wave func- 
tion that resides strictly within the lowest LL. For elec- 
trons, at least two ways of accomplishing this have been 
useful; 2 -' 3 — they both give lowest LL wave functions that 
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are very close to the exact eigenstates. The more con- 
venient of these two methods^ relies on having at least 
two factors of $1; it does not work for but requires 

$f$n- That is the reason for defining the projection as 
in Eq. I|15fl . We refer the reader to the literature^ 3 , for the 
explicit construction of the lowest Landau level projected 
wave functions Vlll^i^u, which can be used here with- 
out change. The presence of $1 in the denominator is 
not a cause for concern, because Vlll&i^u, being anti- 
symmetric, also contains the factor <£>! in it. We have not 
tested the relative merits of this method of projection as 
opposed to a direct projection, but, based on our experi- 
ence with fermions, we expect them to produce more or 
less the same lowest LL wave function. 

To compare ^ B with the exact wave functions 4 ,e2: , we 
will calculate their overlap: 



O 2 



ex \yX,B \|2 

/(„+l)l^„/( n+ l)/l 



I (K X /(n + l) I K X /(n+l))\ I <*n/(n+l)l*n/(n+l))l ' 



(16) 

For the Metropolis Monte Carlo evaluation, it is conve- 
nient to rewrite it as 



O 



2 _ 



l^/(nH 



>I/K*£ 



/(n+l) 



n/(n+l) 



*n/(n+l) 



B 

n/(n+l) 



*n/(n+l)/ 



Then, using the wave function ty B 



/(r. 



(17) 

^ as the sampling 



weight, both the numerator and the denominator can 
be calculated simultaneously. ^n/(n+i) represents either 
the ground state wave function or the CF exciton wave 
function at v = n/(n+l). The corresponding exact wave 
functions are obtained by the Lanczos algorithm. 

Another measure of the quantitative accuracy of the 
CF description is the comparison between the predicted 
energy with the exact one. The CF prediction for the 
energy of the ground or excited states is given by 



E = 



,/(r, 



-l)\v\* 



n/(n+l)/ 



\ n/(n+l)l n/(n+l)/ 



(18) 



Even though the wave functions are rather complicated, 
the integral can be evaluated by the Metropolis Monte 
Carlo method. We find it convenient to write the numer- 
ator as 



E^«/(n+l)l 5(2) ^- n i)l*n/(n+l) 



N(N — 1) 



N 

j n^i^ (2) («i-«2)i^ /( „ +1) (ni,n2,. 



1 N(N — 1) 

2 AttR 2 



1=1 

N 



where, in the spherical geometry, the unit vector S~2^ = 
(sin 6i cos ipi , sin 9i sin (pi , cos Oi ) describes the position of 
particles on the surface of the sphere, and we have used 



J df2i/47ri? 2 = 1 in the last step (R = the radius of 
the sphere), which expresses the integral in a form where 
|\I>n/(n+i)({^]-)l 2 can be used as the sampling function. 

For the Monte Carlo evaluation of the overlap, occupa- 
tion number basis states are transformed into real space 
basis wave functions, which are permanents. The perma- 
nent is an analog of a determinant of a square matrix M 
with elements Mjk in which all signs are taken as positive 
in the expansion of minors. In general, it can be written 
as per(M J 7 C ) = Tlj=i -^KjXj)' w here a are permuta- 
tions of N indices. We evaluate the permanents using 
the Ryser algorithm. 39 Typically we perform 10 3 ~ 10 D 
iterations in one Monte Carlo run. For larger systems, 
the majority of the computational time is spent on eval- 
uating permanents. For example, we need to evaluate 
61108 permanents at each Monte Carlo step for a system 
of N = 12 bosons at v = 3/4 which takes approximately 
480 CPU hours for 10 3 iterations on a single node of a PC 
cluster, consisting of dual 2.4GHz PentiumlV processors, 
to accumulate the desired accuracy. We use as many as 
10 nodes to increase the efficiency. 

In the energy calculation, the wave function in Eq. I|15|) 
consists of a linear combination of several determinants 
for an excited state at a given angular momentum L. 
(For the ground state, only one determinant needs to be 
evaluated.) The calculation for energy is far less time 
consuming than that for the overlap. We perform about 
1.2 x 10 7 iterations in a single Monte Carlo run. The 
quoted statistical uncertainty in the calculation reflects 
one standard deviation from 10 independent runs. To 
give an idea of the computation time, approximately 40 
CPU hours are needed for the ground state energy of v = 
3/4 at N = 12. In Eq. i fHty . the positions of the first two 
particles are identical. To avoid numerical division by 
zero, we set fii = SI2 + Sfl. The results are independent 
of \5fl\ provided it is sufficiently small; we typically use 
|<5fl| = 10- 6 . 

The Moore- Read wave function is known to be the ex- 
act ground state for a three body interactio n 25 ! 40 



Hpfaff 



E * (2 

i<j<k 



(20) 



It can therefore be obtained by exact diagonalization us- 
ing Lanczos algorithm in the spherical geometry. That 
provides a direct method to evaluate the scalar product 
involved in the overlap calculation. 



V. QUANTITATIVE COMPARISONS 



(19) 



The results of our study are summarized in Table I and 
Figs. 1-3. The Table I gives the overlaps of exact eigen- 
states ^ ex with the wave functions ^ B for the ground 
state and first excited state. We make the following re- 
marks. 
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FIG. 1: The low-energy spectrum of (a) TV = 9 and (b) 
N = 10 interacting bosons at v — 1/2, interacting with a 
delta function interaction with strength g. Dashes represent 
exact results, while the dots show predictions of the composite 
fermion theory. Spherical geometry is used in the calculation, 
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FIG. 2: The low-energy spectrum of (a) N = 10 and (b) 
N = 12 interacting bosons at v = 2/3. Various symbols have 
the same meanings as in Fig. 1. 



(i) For v = 1/2, 2/3, and 3/4, the structure of the 
low-energy states is in clear correspondence with that 
of fermions at v* = 1, 2, and 3. In all these cases, 
the ground state is a uniform state (L — 0), well sepa- 
rated from the other states by a gap, as shown earlier iSS 
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FIG. 3: The low-energy spectrum of (a) TV = 9 and (b) TV = 
12 interacting bosons at v — 3/4. See the caption of Fig. 1 
for the definition of various symbols. 



The quantum numbers of the low-energy excitations at 
v = 1/2 and 2/3 can also be understood by analogy to 
v* = 1 and 2, although the correspondence for excita- 
tions is poor for v = 3/4. Thus, an inspection of the 
low-energy spectrum of interacting bosons at v clearly 
shows similarity with fermions at v* . 

(ii) At v = 1/2, the wave function <I>f is known to be 
exact. Our calculations of the overlap and energy for this 
filling constitute a non-trivial test of the correctness of 
our computer codes. 

(iii) The low-energy excited states at v = 1/2 are ex- 
tremely well described, quantitatively, as excitons of com- 
posite fermions The composite- fermion theory pre- 
dicts that there is a single state at orbital angular mo- 
menta from L = 2 to L = TV, which is clearly the case 
in exact diagonalization studies. (At u* = 1 there is also 
an exciton at L = 1, but its wave function is annihilated 
upon projection into the lowest Landau level^i) Further, 
the energy of the CF exciton is in excellent agreement 
with the exact energy. Previous studie s 20 ' 29 studied the 
excited states at v — 1/2 by exact diagonalization, but 
did not provide a microscopic understanding. 

(iv) For v = 2/3, the CF theory provides an excellent 
approximation for the ground state, with very high over- 
laps and very good energies for 10 and 12 particles. The 
CF theory again correctly predicts the quantum numbers 
of the low energy excitations, and also the qualitative 
shape of the exciton dispersion curve, but the energies 
are now off by up to ~ 50%. At v = 3/4, the situa- 
tion becomes worse. In accordance with the prediction 
of the CF theory, the ground state has L — 0, but no 
well defined branch of excitations may be identified with 
the CF exciton branch; furthermore, the energies pre- 
dicted by the CF theory are quite inaccurate, for both 
the ground and excited states. These studies show that 
the CF description worsens with increasing n along the 
sequence v = n/(n + 1). 

(v) At v = 1, a good account of the ground state is 
obtained through analogy to a paired fermion state, as 
can be seen from the overlaps given in the last column 
of the Table I. This result is consistent with the earlier 
studies in the toroidal and disk geometriesi&i£*i&. 

One may ask to what extent the difference between 
bosons at v = n/(n+ 1) and electrons at v e — nj (2n + 1) 
has to do with the fact that the bosons are interacting 
via a short-range, contact interaction, as opposed to the 
long-range Coulomb interaction for the electrons. To in- 
vestigate this issue, we obtain exact wave functions for 
a system of charged bosons interacting via the Coulomb 
interaction. Table II presents their overlaps with vari- 
ous wave functions. The CF theory is in better agree- 
ment with the Coulomb states at v — nj(n + 1), but 
the overall behavior is qualitatively similar. For exam- 
ple, the overlaps for v — 3/4 are not high, and much 
smaller than those at v e = 3/7 for the electron FQHE. 
The paired wave function is also a better approximation 
for the Coulomb ground state at v = 1 than it is for 
the hard-core interaction; in contrast, it is not valid for 
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TABLE I: The overlap of the exact wave functions for the ground state and the first excited state at v — 1/2, 2/3, 3/4, and 
1 with the trial wave functions of Laughlin [y — 1/2 ground state), Moore and Read (y — 1 ground state), and Jain (other 
states), for several particle numbers N. O gr is the overlap for the ground state, and O ex for the first excited state, which occurs 
at the orbital angular momentum L. The definition of the overlap is given in Eq. 11611 . The statistical uncertainty in the last 
two digits is shown in parentheses when it is larger than 10~ . For v = 1, only the ground state overlap is shown, which has 
been evaluated exactly. 
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TABLE II: The overlap of the Coulomb ground state wave function at v — 1/2, 2/3, 3/4, and 1 with the trial wave functions of 
Laughlin (y = 1/2 ground state), Moore and Read (v — 1 ground state), and Jain [y = 2/3 and 3/4 ground states), for several 
particle numbers N. (Table I dealt with the exact wave functions for a short range interaction.) The statistical uncertainty in 
the last two digits is shown when it is larger than 10 -5 . The overlaps for v = 1 (last column) are exact. 



the electronic state at the corresponding filling v e = 1/2. 
These observations indicate that both the particle statis- 
tics and the nature of the interaction are responsible for 
the differences in the behaviors of fermions and bosons 
in a magnetic field. 



VI. CONCLUSION 

The above results allow us to make the following con- 
clusions: (i) The mapping into free fermions is quali- 
tatively valid for a range of parameters. It correctly 
captures the incompressibility of the ground state at 
v = 1/2, 2/3 and 3/4. (ii) The mapping is also quantita- 
tively very accurate for the ground state and excitations 
at v = 1/2 and for the ground state at v = 2/3, but be- 
comes progressively worse with increasing n. This implies 
that the residual interactions between fermions become 
increasingly more important with n, and must be con- 
sidered for a more complete understanding, (iii) A qual- 
itative indication of the breakdown of the /ree-fermion 
model is the appearance of a paired state at v = 1/2; 
noranteracting fermions would have produced a Fermi 
sea here. The residual interactions between fermions thus 
cause a qualitative change in the nature of the state be- 
yond certain n, presumably through a pairing instability; 



we cannot, however, ascertain from our study for what n 
a phase transition occurs. 

Overall, these results establish that the mapping into 
strictly free fermions is valid only for a limited range 
of parameters, but the mapping into weakly interact- 
ing fermions provides a useful starting point for a wider 
range of parameters. One might ask why the interacting 
fcrmion language is to be preferred over the original in- 
teracting boson model. The reason is that the interaction 
between the fermions is much weaker, with a large part 
of the repulsive interaction taken care of by the Pauli 
avoidance. 

The bosonic FQHE should be contrasted with the 
FQHE of electrons along the sequence v e = n/(2n + 1), 
for which the mapping into free composite fermions re- 
mains valid for the entire parameter range. The wave 
functions of Eq. © are known to provides an excellent 
description of the state for 1/3, 2/5, and 3/7, where exact 
results are available and presumably also for higher 
n, as evidenced by the experimental observation of many 
fractions along the sequence as well as of the Fermi sea 
of composite fermions at v = 1/2«22j2! 

It is noteworthy that pairing occurs in a model of 
bosons with purely repulsive interactions. As stressed 
in Refl38t even when the interaction between the original 
particles is repulsive, the effective interaction between the 
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emergent particles may be attractive. That appears to 
be the case at v = 1. Here, bosons dress themselves with 
vortices to turn into fermions, but that presumably over- 
compensates for the repulsive interaction, thereby pro- 
ducing an attractive interaction between the fermions. 
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